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Abstract 

We consider the system of Maxwell-Stefan equations which describe mul- 
f} , ticomponent diffusive fluxes in non-dilute solutions or gas mixtures. We 



apply the Perron-Frobenius theorem to the irreducible and quasi-positive 
matrix which governs the flux-force relations and are able to show normal 
<**] . ellipticity of the associated multicomponent diffusion operator. This pro- 

vides local-in-time wellposedness of the Maxwell-Stefan multicomponent 
diffusion system in the isobaric, isothermal case. 
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r^- ■ 1 Introduction 

O " 

On the macroscopic level of continuum mechanical modeling, fluxes of chemical 
components (species) are due to convection and molecular fluxes, where the 
latter essentially refers to diffusive transport. The almost exclusively employed 
constitutive "law" to model diffusive fluxes within continuum mechanical models 
Vh ' is Fick's law, stating that the flux of a chemical component is proportional to the 

gradient of the concentration of this species, directed against the gradient. There 
is no influence of the other components, i.e. cross-effects are ignored although 
well-known to appear in reality. Actually, such cross-effects can completely 
divert the diffusive fluxes, leading to so-called reverse diffusion (up-hill diffusion 
in direction of the gradient) or osmotic diffusion (diffusion without a gradient). 
This has been proven in several experiments, e.g. in a classical setting by Duncan 
and Toor; see [7]. 

To account for such important phenomena, a multicomponent diffusion ap- 
proach is required for realistic models. The standard approach within the the- 
ory of Irreversible Thermodynamics replaces Fickian fluxes by linear combi- 
nations of the gradients of all involved concentrations, respectively chemical 



potentials. This requires the knowledge of a full matrix of binary diffusion coef- 
ficients and this diffusivity matrix has to fulfill certain requirements like positive 
semi-definiteness in order to be consistent with the fundamental laws from ther- 
modynamics. The Maxwell-Stefan approach to multicomponent diffusion leads 
to a concrete form of the diffusivity matrix and is based on molecular force 
balances to relate all individual species velocities. While the Maxwell-Stefan 
equations are successfully used in engineering applications, they seem much less 
known in the mathematical literature. In fact we are not aware of a rigorous 
mathematical analysis of the Maxwell-Stefan approach to multicomponent diffu- 
sion, except for [5] which mainly addresses questions of modeling and numerical 
computations, but also contains some analytical results which are closely related 
to the present considerations. 

2 Continuum Mechanical Modeling 
of Multicomponent Fluids 

We consider a multicomponent fluid composed of n chemical components Aj. 
Starting point of the Maxwell-Stefan equations are the individual mass balances, 
i.e. 

d tPi + 6Jy{p i u i ) = Rf\ (1) 

where pt = Pi(t,y) denotes the mass density and u^ = Uj(t, y) the individual 
velocity of species Aj. Note that the spatial variable is denoted as y, while the 
usual symbol x will refer to the composition of the mixture. The right-hand side 
is the total rate of change of species mass due to all chemical transformations. 
We assume conservation of the total mass, i.e. the production terms satisfy 
S"=i ^i 0t = 0. Let p denote the total mass density and u be the bary centric 
(i.e., mass averaged) velocity, determined by 

n n 

P~^2pu pu:=^2 Pl u t . 

i=l z=l 

Summation of the individual mass balances ([T]) then yields 

d t p + div (pu) = 0, (2) 

i.e. the usual continuity equation. 

In principle, a full set of n individual momentum balances should now be 
added to the model; cf. [11 . But in almost all engineering models, a single set 
of Navier-Stokes equations is used to describe the evolution of the velocity field, 
usually without accounting for individual contributions to the stress tensor. One 
main reason is a lack of information about appropriate constitutive equations 
for the stress in multicomponent mixtures; but cf. [16] . For the multicomponent, 
single momentum model the barycentric velocity u is assumed to be determined 
by the Navier-Stokes equations. Introducing the mass diffusion fluxes 

ji :=pi(ui-u) (3) 

and the mass fractions Yi := pi/ p, the mass balances ([T]) can be rewritten as 

pdtY, + pu -VY t + div h = Rl ot . (4) 



In the present paper, main emphasis is on the aspect of multicomponent dif- 
fusion, including the cross-diffusion effects. Therefore, we focus on the special 
case of isobaric, isothermal diffusion. The (thermodynamic) pressure p is the 
sum of partial pressures pi and the latter correspond to CiRT in the general 
case with Ci denoting the molar concentration, R the universal gas constant and 
T the absolute temperature; here Cj = pi/Mi with Mi the molar mass of species 
Ai. Hence isobaric conditions correspond to the case of constant total molar 
concentration ctot, where ctot : = X)"=i c *- Still, species diffusion can lead to 
transport of momentum because the Mj are different. Instead of u we therefore 
employ the molar averaged velocity defined by 

n 

CtotV := y^CjUi- (5) 

Note that other velocities are used as well; only the diffusive fluxes have to be 
adapted; see, e.g., [2D]. With the molar averaged velocity, the species equations 
([T]) become 

dta + div (c 4 v + J,) = r| ot (6) 

with r* ot := i?* ot /Mi and the diffusive molar fluxes 

J. t := c, : (u, - v). (7) 

Below we exploit the important fact that 

n 

E J * = °- ( 8 ) 

x=l 

As explained above we may now assume v = in the isobaric case. In this 
case the species equations (JDJ) simplify to a system of reaction-diffusion systems 
given by 

d tCi + div 3i = rt°\ (9) 

where the individual fluxes 3i need to be modeled by appropriate constitutive 
equations. The most common constitutive equation is Fick's law which states 
that 

Ji = -Degrade* (10) 

with diffusivities Di > 0. The diffusivities are usually assumed to be constant, 
while they indeed depend in particular on the composition of the system, i.e. 
Di = A(c) with c := (ci . . . , c n ). Even if the dependence of the Di is taken 
into account, the above definition of the fluxes misses the cross-effects between 
the diffusing species. In case of concentrated systems more realistic constitutive 
equations are hence required which especially account for such mutual influences. 
Here a common approach is the general constitutive law 



j'=i 



A j grade,- (11) 



with binary diffusivities D^ = Dy(c). Due to the structure of the driving forces, 
as discussed below, the matrix D = [Dij] is of the form D(c) = L(c)G"(c) 
with a positive definite matrix G"(c), the Hessian of the Gibbs free energy. 
Then, from general principles of the theory of Irreversible Thermodynamics, it 
is assumed that the matrix of transport coefficients L = [Lij] satisfies 

• L is symmetric (the Onsager reciprocal relations) 

• L is positive semidefinite (the second law of thermodynamics). 
Under this assumption the quasilinear reaction-diffusion system 

<9 t c + div (-D(c) Vc) =r(c), (12) 

satisfies - probably after a reduction to n—1 species - parabolicity conditions suf- 
ficient for local-in-time wellposedness. Herer(c) is short for (r\ ot (c), . . . ,r^ ot (c)). 

A main problem now is how realistic diffusivity matrices together with their 
dependence on the composition vector c can be obtained. 

Let us note in passing that Herbert Amann has often been advocating that 
general flux vectors should be considered, accounting both for concentration 
dependent diffusivities and for cross-diffusion effects. For a sample of his con- 
tributions to the theory of reaction-diffusion systems with general flux vectors 
see pQ, [2] and the references given there. 

3 The Maxwell- Stefan Equations 

The Maxwell-Stefan equations rely on inter-species force balances. More pre- 
cisely, it is assumed that the thermodynamical driving force d^ of species Ai 
is in local equilibrium with the total friction force. Here and below it is often 
convenient to work with the molar fractions Xi :— c^/ctot instead of the chemical 
concentrations. From chemical thermodynamics it follows that for multicompo- 
nent systems which are locally close to thermodynamical equilibrium (see, e.g., 
|20j ) the driving forces under isothermal conditions are given as 

d ' = EY gmdAil ^ 

with /ii the chemical potential of species Ai . Equation (|13[) requires some more 
explanation. Recall first that the chemical potential /Zj for species Ai is defined 

as 

where G denotes the (volume-specific) density of the Gibbs free energy. The 
chemical potential depends on c, , but also on all other Cj as well as on pressure 
and temperature. In the engineering literature, from the chemical potential a 
part /i® depending on pressure and temperature is often separated and, depend- 
ing on the context, a gradient may be applied only to the remainder. To avoid 
confusion, the common notation in use therefore is 

dm 9/Xj 



Here Vr^fJ-i means the gradient taken under constant pressure and tempera- 
ture. In the isobaric, isothermal case this evidently makes no difference. Let us 
also note that G is assumed to be a convex function of the Ci for single phase sys- 
tems, since this guarantees thermodynamic stability, i.e. no spontaneous phase 
separations. For concrete mixtures, the chemical potential is often assumed to 
be given by 

Hi = /i-+ RT In a t (15) 

with Oj the so-called activity of the z-th species; equation (|15[) actually implicitly 
defines a». In (|15p . the term $ depends on pressure and temperature. For a 
mixture of ideal gases, the activity a, equals the molar fraction cc,. The same 
holds for solutions in the limit of an ideally dilute component, i.e. for Xi — > 0+. 
This is no longer true for non-ideal systems in which case the activity is written 
as 

a l = -f i Xi (16) 

with an activity coefficient 7$ which itself depends in particular on the full 
composition vector x. 

The mutual friction force between species i and j is assumed to be propor- 
tional to the relative velocity as well as to the amount of molar mass. Together 
with the assumption of balance of forces this leads to the relation 

d; = - Y^ hi x i x i ( u * - u 3 ) ( 1 7 ) 

with certain drag coefficients fy > 0; here / y = fji is a natural mechanical 
assumption. Insertion of (fT5]) and introduction of the so-called Maxwell- Stefan 
(MS) diffusivities By = 1/fij yields the system 

— ^grad/i, = - V J - ' J for i = 1, . . . , n. (18) 

RT *-■ '. ctotBy 

3=£* 

The set of equations (|T5)) together with (JSJ) forms the Maxwell- Stefan equations 
of multicomponent diffusion. The matrix [By] of MS-diffusivities is assumed to 
be symmetric in accordance with the symmetry of [fij]. Let us note that for 
ideal gases the symmetry can be obtained from the kinetic theory of gases; cf. 
[9] and [14] . The MS-diffusivities D^ will in general depend on the composition 
of the system. 

Due to the symmetry of [D,j], the model is in fact consistent with the On- 
sager reciprocal relations (cf. [18J as well as below), but notice that the D^ are 
not to be inserted into (fTTj) . i.e. they do not directly correspond to the D y there. 
Instead, the MS equations have to be inverted in order to provide the fluxes J^ . 

Note also that the Ansatz (fTT)) implies ^\ d^ = because of the symmetry 
of [fij], resp. of [By]. Hence Y^i^i = is necessary in order for (TT71) to be 
consistent. It in fact holds because of (and is nothing but) the Gibbs-Duhem 
relation, see e.g. [12] ■ The relation J^ i d^ = will be important below. 

Example (Binary systems). For a system with two components we have 
di(=-d 2 ) = (jca^-a^jjj). (19) 

CtotBi2 



Using x\ + X2 — 1 and Ji + J2 = one obtains 
Ji(= -J2) = --^cigrad/xi. 



(20) 



Writing c and J instead of c\ and Ji, respectively, and assuming that the chem- 
ical potential is of the form [i — fjP + RT\a{^c) with the activity coefficient 
7 = 7(c) this finally yields 



J = -D 



12 (l + ^)gradc. 

V -VIC / 



(21) 



Inserting this into the species equation leads to a nonlinear diffusion equation, 
namely 



d t c- A<t>(c) = r(c), 



(22) 



where the function <fi : R — > R satisfies <f>'(s) = Di2(l + S7'(s)/7(s)) and, 
say, </>(0) = 0. Equation (f2"2"|) is also known as the filtration equation (or, the 
generalized porous medium equation) in other applications. Note that well- 
known pde-theory applies to (|22[) and especially provides well-posedness as soon 
as 4> is continuous and nondecreasing; cf., e.g., [21] . The latter holds if s — > sj(s) 
is increasing which is nothing but the fact that the chemical potential fi of 
a component should be an increasing function of its concentration. This is 
physically reasonable in systems without phase separation. 



4 Inversion of the Flux- Force Relations 

In order to get constitutive equations for the fluxes J^ from the Maxwell-Stefan 
equations, which need to be inserted into (O, we have to invert (JT5J). Now 
(fT8| alone is not invertible for the fluxes, since these are linearly dependent. 
Elimination of J„ by means of © leads to the reduced system 



Ctot 



di 



B 



Ji 



(23) 



where the (n — 1) x (n — l)-matrix B is given by 



Bij - { 



n77 + 2 d^T for * = j ^ with Xn = l ~ ^ m< ™ Xm ^' 



(24) 



Assuming for the moment the invertibility of B and letting fii be functions of 
the composition expressed by the molar fractions x = (xi, . . . ,x„), the fluxes 
are given by 



Ji 



V«i 



= -c to tB" 



(25) 



Vav 



where 

din 7, 
dxj 

captures the thermodynamical deviations from the ideally diluted situation; here 
8ij denotes the Kronecker symbol. 



T = [Tij] with Tij = Sij + Xi -^ (26) 



Example (Ternary systems). We have 



B 



i), ; ; C2 ln^ d 13 J Xi \b 12 d 1:; 



(27) 



and dct(B - tl) = t 2 - tr B t + dot B with 

detB = h - 1 - — >min{ , — - , } (28) 

D12D13 D12D23 D13D23 D12D13 D12D23 D13D23 

and 

, „ Xi + x 2 Xi + x 3 x 2 + x 3 . 1 1 1 

trB = 1 1 — > 2mm{ , , }. (29) 

012 £>13 E*23 t>12 013 023 

It is easy to check that (tr B) 2 > 3 det B for this particular matrix and there- 
fore the spectrum of B _1 is in the right complex half-plane within a sector of 
angle less than tt/6. This implies normal ellipticity of the differential operator 
B _1 (x)(— Ax). Recall that a second order differential operator with matrix- 
valued coefficients is said to be normally elliptic if the symbol of the principal 
part has it's spectrum inside the open right half-plane of the complex plane; see 
section 4 in |2J for more details. This notion has been introduced by Herbert 
Amann in [I] as the appropriate concept for generalizations to more general 
situations with operator-valued coefficients. 

Consequently, the Maxwell-Stefan equations for a ternary system are locally- 
in-time wellposed if T = I, i.e. in the special case of ideal solutions. The latter 
refers to the case when the chemical potentials are of the form (fl~5|) with 7; = 1 
for all i. Of course this extends to any T which is a small perturbations of I, 
i.e. to slightly non-ideal solutions. 

Let us note that Theorem 1 below yields the local-in-time wellposedness also 
for general non-ideal solutions provided the Gibbs energy is strongly convex. 
Note also that the reduction to n — 1 species is the common approach in the 
engineering literature, but invertibility of B is not rigorously checked. For n = 4, 
the 3 x 3-matrix B can still be shown to be invertible for any composition due 
to Xi > and YljXj = 1. Normal ellipticity can no longer be seen so easily. 
For general n this approach is not feasible and the invariant approach below is 
preferable. 

Valuable references for the Maxwell-Stefan equations and there applications 
in the Engineering Sciences are in particular the books [1], [5], [20] and the 
review article |12) . 

5 Wellposedness of the Maxwell- Stefan equations 

We first invert the Maxwell- Stefan equations using an invariant formulation. 
For this purpose, recall that ^\ u, = holds for both m = Ji and Ui = d,. We 



therefore have to solve 

A3 = c to td in E = {ueR n 

where A — -A(x) is given by 

-Si di 



£■ 



= 0}, 



(30) 



A = 






with Si = 



EXk 

k=ii lL 



A - — 

lJ - D 



'j 



The matrix A has the following properties, where x 3> means Xi > for all i: 

(i) iV(A) = spanjx} for x = (x\, . . . , X n ). 

(ii) i?(A) = {e}^fore = (l,...,l). 

(iii) A = [aij] is quasi-positive, i.e. a,ij > for i =/= j. 

(iv) If x S> then A is irreducible, i.e. for every disjoint partition I U J of 
{1, . . . , n} there is some (z, j) E I x J such that Ojj 7^ 0. 

Due to (i) and (ii) above, the Perron-Frobenius theorem in the version for quasi- 
positive matrices applies; cf. [TU] or [TTj. This yields the following properties 
of the spectrum cr(A): The spectral bound s(A) := max{ReA : A G &(A)} 
is an eigenvalue of A, it is in fact a simple eigenvalue with a strictly positive 
eigenvector. All other eigenvalues do not have positive eigenvectors or positive 
generalized eigenvectors. Moreover, 

ReA<s(A) forallAG(7(i),A^s(4 

From now on we assume that in the present case x is strictly positive. Then, 
since x is an eigenvector to the eigenvalue 0, it follows that 

a{A) C {0} U{z€d:Rcz< 0}. 

Unique solvability of (|3"0)) already follows at this point. In addition, the same 
arguments applied to Ay, :— A — /i(x ® e) for /1 G R yield 

cr{Ay) C {-A*} U {z Gd : Re z < -ju} for all small fi > 0. 

In particular, A^ is invertible for sufficiently small /1 > and 

J = -c tot (A-/i(x®e)) _1 d (31) 

is the unique solution of (|3"0")) . Note that A^y = d with die implies yle 
and Ay = d. A similar representation of the inverted Maxwell-Stefan equations 
can be found in [5]. 

The information on the spectrum of A can be significantly improved by 
symmetrization. For this purpose let X = diag(xi, . . . x n ) which is regular due 
to x > 0. Then A s := X^ Axi satisfies 



As 



Si 



'1.1 
-s n 



D ifc ' 



k^i 






i.e. As is symmetric with N(As) — spanj-y/x}, where y/x i := ^fxi. Hence the 
spectrum of As and, hence, that of A is real. Moreover, 

As(a) = A s — ocy/x ® V* 

has the same properties as As for sufficiently small a > 0. In particular, A5 is 
quasi-positive, irreducible and -y/x 3> is an eigenvector for the eigenvalue —a. 
This holds for all a < 5 := min{l/Dy : i ^ j}. Hence we obtain the improved 
inclusion 

(7(A) \ {0} = a(A s (a)) \ {-a} for all a £ [0,6). 

Therefore 

a {A) C (-00, -<5] U {0}, (32) 

which provides a uniform spectral gap for A sufficient to obtain normal ellipticity 
of the associated differential operator. 

In order to work in a subspace of the composition space R™ instead of a 
hyperplane, let m — Ci — c® ot /n such that J^ q = const is the same as u G E = 
{u £ R ra : Y^i u i — 0}- Above we have shown in particular that A\ E : E — > E is 
invertible and 

[J,] = X*^)- 1 *-* [d,] - ±,xi(A S[i: )-'xi [V M i] (33) 

with the symmetrized form As of A and E := X?E = {y/x}- 1 . Note that this 
also shows the consistency with the Onsager relations. To proceed, we employ 
(TTJI to obtain the representation 

[J^xHA^r'xiG"^)^. (34) 

Inserting (|34[) into @ and using c tot Xi = Ui + c® ot /n, we obtain the system of 
species equations with multicomponent diffusion modeled by the Maxwell-Stefan 
equations. Without chemical reactions and in an isolated domain Q C R™ (with 
v the outer normal) we obtain the initial boundary value problem 

d t u + div (-T)(u)Vu) = 0, d v u\ 9 a = 0, u\ t= o = u , (35) 

which we will consider in D>(Q.; E). Note that Xi{A sl& )- l XiG"(x) from ([541) 
corresponds to — D(u) here. 

Applying well-known results for quasilinear parabolic systems based on i p - 
maximal regularity, e.g. from [3] or |15j . we obtain the following result on local- 
in-time wellposedness of the Maxwell- Stefan equations in the isobaric, isother- 
mal case. Below we call G € C 2 (V) strongly convex if G"(x) is positive definite 
for all xeV. 

Theorem 1 Let O C R with N > 1 be open bounded with smooth dfl. Let 

2— - 

p > ^ and Mo € W p p (f2; E) such that c- > in !1 and c® ot is constant in 
f2. Let the diffusion matrix D(u) be given according to (|34[) . i.e. by 

D(u) = X2(A s| ^) _1 X^G"(x) OTt/i c tot a; 4 = u, + cj? ot /n, 

where G : (0, 00)" — > R is smooth and strongly convex. Then there exists - locally 
in time - a unique strong solution (in the L p -sense) of (I35[) . This solution is in 
fact classical. 



Concerning the proof let us just mention that 

div (-D(ii)Vu) = D(u) (—Am) + lower order terms, 

hence the system of Maxwell-Stefan equations is locally-in-time wellposed if 
the principal part D(u) (—Am) is normally elliptic for all u £ E such that 
c(u) :— u + Ct ot e is close to c°. The latter holds if, for some angle 9 £ (0, ?), 
the spectrum of D(u) £ C{E) satisfies 

ct(D(u)) C S e := {A e(T \ {0} : |arg A| < 9} (36) 

for all u £ E such that \c(u) - c ^ < e for e := miniC°/2, say. For such an 
u £ E, let A ed and v £ E be such that D(w) v — Xv. Let x :— c(u) / Ctot(u) £ 
(0, oo) n and X — diag(xi, . . . , x n ). Then 

X5(A s| ^)- 1 liG"(x)t) = Xv. 

Taking the inner product with G"(x) v yields 

((A sl z)- 1 X 1 zG"( X )v,X 1 zG"( X )v) = X(v,G"(x)v). 

Note that XzG"(x)v £ {^/x}^, hence the left-hand side is strictly positive 
due to the analysis given above. Moreover (v, G"(x) v) > since G is strongly 
convex, hence A > 0. This implies (|36|) for any 9 £ (0, ^) and, hence, local- in- 
time existence follows. 



6 Final Remarks 

A straight-forward extension of Theorem 1 to the inhomogeneous case with 
locally Lipschitz continuous right-hand side / : R™ — » R", say, is possible if 
f(u) £ E holds for all u. Translated back to the original variables (keeping the 
symbol /) this yields a local-in-time solution of 

d t c + div(-D(c)Vc) =/(c), d v c\ da = 0, C| i=0 = c 

for appropriate initial values Co . Then a natural question is whether the solution 
stays componentwise nonnegative. This can only hold if / satisfies 

fi(c) > whenever c > with c, = 0, 

which is called quasi-positivity as in the linear case. In fact, under the considered 
assumption, quasi-positivity of / forces any classical solution to stay nonnegative 
as long as it exists. The key point here is the structure of the Maxwell-Stefan 
equations (|T51) which yields 

J. ( = -A(c)gradci + a F,(c, grade) 

with 

D i (c)=l/J2^r and F,(c,gradc)=A(c)^— J,-. 



10 



Note that Di(c) > and J; becomes proportional to gradc^ at points where a 
vanishes, i.e. the diffusive cross-effects disappear. Moreover, it is easy to check 
that 

div Jj = Di(c)Aci > if q = and grade; = 0. 

To indicate a rigorous proof for the nonnegativity of solutions, consider the 
modified system 

dta + div J, ; (c) = fi(t,c + ) + e, d v C\ d n = 0, C| t=0 = c + ee, (37) 

where r + := max{r, 0} denotes the positive part. Assume that the right-hand 
side / is quasi-positive and that (j3"7) has a classical solution c e for all small e > 
on a common time interval [0,T). Now suppose that, for some i, the function 
rrii(t) — min g Q c|(t,y) has a first zero at to € (0,T). Let the minimum of 
c|(ioj be attained at yo and assume first that yo is an interior point. Then 
cl(t ,yo) = 0, 8scf(t ,yo) > 0, gradc|(i ,yo) = and Ac|(t ,yo) > yields a 
contradiction since fi(to,cf(to,yo)) > 0. Here, because of the specific boundary 
condition and the fact that tt has a smooth boundary, the same argument works 
also if yo is a boundary point. In the limit e — > 0+ we obtain a nonnegative 
solution for e = 0, hence a nonnegative solution of the original problem. This 
finishes the proof since strong solutions are unique. 

Note that non-negativity of the concentrations directly implies L°°-bounds 
in the considered isobaric case due to < Ci < Ctot = c^ ot , which is an important 
first step for global existence. 

The considerations in Section [S] are helpful to verify that the Maxwell-Stefan 
multicomponent diffusion is consistent with the second law from thermodynam- 
ics. Indeed, (|33l) directly yields 

- [J,] : [V W ] = ^(V^r 1 ^ I V ^] ) : (x* [Vim] ) > o, 

i.e. the entropy inequality is satisfied. The latter is already well-known in the 
engineering literature, but with a different representation of the dissipative term 
using the individual velocities; cf. |18j . 

For sufficiently regular solutions and under appropriate boundary conditions 
the entropy inequality can be used as follows. Let V(x) — J n G(x) cfac with G 
the Gibbs free energy density. Let 

W(x, Vx) = - J [Ji] : [Vfu] dx > 0. 
Jn 

Then (V, W) is a Lyapunov couple, i.e. 

V(x(f)) + / W(x(s) 7 \7x(s))ds < V(x(0)) for t > 



and all sufficiently regular solutions. For ideal systems this yields a priori bounds 
on the quantities |Vci| 2 /ci, hence, equivalently, i^-bounds on V-i/cJ. This type 
of a priori estimates is well-known in the theory of reaction-diffusion systems 
without cross-diffusion; see [S], [5] and the references given there for more details. 

In the present paper we considered the isobaric and isothermal case because it 

11 



allows to neglect convective transport and, hence, provides a good starting point. 
The general case of a multicomponent flow is much more complicated, even in 
the isothermal case. This case leads to a Navier-Stokes-Maxwell-Stefan system 
which will be studied in future work. 

Acknowledgement. The author would like to express his thanks to Jan Priiss 
(Halle- Wittenberg) for helpful discussions. 
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